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SECTION  I 


INTRODUCTION 

The  sublimation  of  carbon  may  be  divided  into  several  regimes  depend- 
ing upon  whether  the  solid-vapor  and  gas  phase  reactions  are  in  equilibrium. 
Nonequilibrium  effects  for  both  the  gas  phase  reactions  and  the  solid  vapori- 
zation reaction  are  significant  only  at  low  pressures.  Since  the  intended 
application  of  this  work  is  for  high  pressure  planetary  entry  calculations, 
only  the  equilibrium  case  is  considered.  That  is,  all  gas  phase  reactions 
are  presumed  to  be  in  equilibrium,  and  the  solid  carbon  surface  is  assumed 
to  be  in  equilibrium  with  each  of  the  carbon  sublimation  species. 

The  equilibrium  sublimation  of  carbon  has  been  extensively  investigated 
in  the  literature  for  various  assumptions  concerning  the  sublimation  species 
and  the  dominant  gas  phase  reactions.  In  particular,  the  use  of  differing 
vapor  pressure  laws  (for  the  sublimation  species)  and  associated  thermo- 
chemical data  (heats  of  formation  and  free  energy  functions  for  each  species) 
leads  to  a considerable  variation  in  predicted  surface  mass  transfer  coeffi- 
cients (B').  The  Bibliography  is  by  no  means  complete  but  is  typical  of  the 
available  work. 

The  one  common  feature  of  the  referenced  work  (with  the  exception  of 
Ref.  1)  is  that  the  surface  mass  transfer  coefficients  were  obtained  by  an 
iterative  numerical  scheme  as  is  usually  employed  in  general  purpose  equi- 
librium chemistry  codes  such  as  ACE  (Ref.  2).  Bartlett  (Ref.  1),  by 

^Bartlett,  E.  P.,  "Analytical  and  Graphical  Prediction  of  Graphite  Ablation 
Rates  and  Surface  Temperatures  During  Reentry  at  25,000  to  45,000  Feet 
per  Second,  " Air  Force  Flight  Test  Center,  Edwards  Air  Force  Base, 
California,  TDR-63-40,  March  1964. 

2 

Powars,  C.A.,  and  R.M.  Kendall,  "User's  ManuaU-Aerotherm  Chemical 
Equilibrium  (ACE)  Computer  Program,  " Aerotherm  Corporation, 

May  1969. 
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assuming  only  one  sublimation  species  (C^)  and  ignoring  carbon-nitrogen 
reactions,  obtained  an  explicit,  analytic  solution  for  the  surface  mass  trans- 
fer coefficient.  For  carbon  sublimation  in  an  atmosphere  of  elemental  car- 
bon, nitrogen,  oxygen,  and  an  inert  (denoted  by  I)  monatomic  species,  it  is 
possible  to  generalize  Bartlett's  work  to  include  any  number  of  carbon 
sublimation  species  and  all  of  the  significant  carbon-nitrogen  reactions, 
while  still  obtaining  an  analytic  solution  for  the  mass  transfer  parameter 
(B').  Similarly,  an  analytic  solution  is  possible  in  an  atmosphere  of  C,  N, 

O,  H and  I. 

Apparently,  it  is  not  generally  (if  at  all)  recognized  that  analytic  solu- 
tions are  possible  for  the  equilibrium  sublimation  of  carbon  in  many  atmos- 
pheres of  interest.  It  is  the  purpose  of  this  report  to  present  these  solutions 
for  the  turface  mass  transfer  coefficient  and  compare  them  with  exact  results 
obtained  from  the  ACE  computer  code. 

Aside  from  their  own  intrinsic  elegance,  there  exist  several  utilitarian 
reasons  for  developing  analytic  solutions  for  B',  First,  precomputed  values 
of  B'  are  valid  only  for  a specified  set  of  species  thermochemical  data  and 
for  an  atmosphere  of  fixed  elemental  composition.  If  the  thermochemical 
data  or  the  atmospheric  composition  is  changed,  the  value  of  B'  must  be 
recalculated.  Analytic  solutions  eliminate  the  need  to  precompute  B',  there- 
by allowing  the  rapid  assessment  of  differing  thermochemical  data  and  atmos- 
pheric composition  upon  the  mass  transfer.  Secondly,  the  solution  of  the 
energy  equation  coupled  to  the  surface  energy  balance  involves  an  iterative 
scheme  which  is  quite  sensitive  to  the  temperature  derivative  of  B'.  This 
derivative  is  diffii'ult  to  evaluate  accurately  from  tabulated  values  of  B', 
but  is  easily  calculated  from  an  analytic  solution. 


I 
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The  analytic  solution  for  carbon  sublimation  in  an  atmosphere  of 
elemental  carbon,  nitrogen,  oxygen,  and  an  inert  species  is  essentially 
exact  and  is  presented  in  Section  II.  The  inclusion  of  hydrogen  into  the 
atmospheric  composition  is  treated  in  Section  III  and  yields  an  analytic 
solution  which  is  slightly  less  precise  than  that  of  Section  II  (due  to  the 
neglect  of  carbon-nitrogen  reactions),  but  which  is  still  in  excellent  agree- 
ment with  exact  solutions. 


( 

I ! 

i; 


i 


i 

I 

i 


i 
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SECTION  II 

CARBON  SUBLI’'4ATION  IN  A C-N-O-I  ATMOSPHERE 

As  discussed  earlier,  the  objective  is  to  obtain  an  analytic  expression 
for  the  surface  mass  transfer  parameter  (B').  The  elemental  composition  of 
the  atmosphere  is  assumed  to  consist  of  carbon,  nitrogen,  oxygen,  and  an 
inert  monatomic  species  (denoted  by  I).  This  specification  includes  a variety 
of  atmospheres  of  interest,  such  as  air  and  a Mars  atmosphere  consisting  of 

N^. 

The  independent  variables  are  the  surface  temperature  and  pressure 
and  the  elemental  composition  at  the  boundary  layer  edge.  By  definition, 
the  elemental  mass  fractions  (denoted  by  a superscript  tilde)  must  sum  to 
unity. 


+ K_  + K,,,  + = 1 

ONI 
e e e 


(1) 


where  the  subscript  e refers  to  boundary  layer  edge  conditions. 

Since  the  solid  carbon  surface  is  assumed  to  be  in  equilibrium  with  its 
vapor,  the  partial  pressure  of  each  sublimation  species  is  equal  to  its  vapor 
pressure.  The  vapor  pressures,  presented  in  Appendix  A,  are  known  func- 
tions of  surface  temperature.  It  is  precisely  this  knowledge  of  the  partial 
pressures  of  the  carbon  species  which  allows  an  analytic  solution.  Thus,  in 
all  that  follows,  the  partial  pressures  of  the  carbon  species  have  been 
replaced  by  their  respective  vapor  pressures. 

A number  of  commonly  employed  assumptions  are  used  to  simplify  the 
problem.  All  the  oxygen  at  the  carbon  surface  is  assumed  to  be  bound  up  in 
CO  and  thus  is  not  available  for  the  formation  of  O^,  NO,  etc.  The 

surface  temperature  is  presumed  to  be  low  enough  so  that  monatomic  oxygen 
and  nitrogen  are  not  present  at  the  wall.  Ionized  species  are  ignored. 


-9- 


r 


i 


An  arbitrary  number  of  carbon  sublimation  species  are  allowed  in  the 
formulation  and  all  of  the  domine'-nt  carbon-nitrogen  reactions  are  included. 
The  carbon-nitrogen  reactions  may  be  expressed  in  any  form  as  long  as  the 
system  is  linearly  independent.  Furthermore,  since  all  the  carbon  sublima- 
tion species  are  in  equilibrium  with  one  another,  the  carbon-nitrogen  reac- 
tions may  be  expressed  in  terms  of  any  one  (or  combination)  of  the  sublima- 
tion species.  With  these  points  in  mind,  the  carbon-nitrogen  reactions  are 
expressed  as 

+ 2C-*2  CN 
C2N2-2  CN 
C^N2-^2  CN  + 2C 

C^N-^CN  + C (2; 

Notice  that  no  dissociation  reactions  involving  the  carbon  sublimation  species 
are  explicitly  included  in  the  formulation.  This  is  not  necessary  because  the 
use  of  the  carbon  vapor  pressure  laws  guarantees  gas  phase  equilibrium 
among  the  sublimation  species. 

The  elemental  mass  balances  at  the  wall  are  given  by 

^ 16 

^Nw  ~ '''  IT  ^CN  52  '''  TS’  ^C^N^  38  ^C2N 

^C  " TS  ^CO  Z€  ^CN  52  ^C^N  tE  ^C  . N 
w 2 2 4 2 

2 1 

Kj  = Kj  (3) 
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where  subscript  w refers  to  conditions  at  the  wall.  The  molecular  weight  is 
given  by 


I 


^CN  , ^2^2  , ^4^2  , ^2^  ^ ^i 

"2^  52  “T?  38  Mj 

i 


Note  that  the  inert  species  is  characterized  solely  by  its  molecular  weight 

(Mj). 

The  equilibrium  relationships  for  the  carbon-nitrogen  reactions  may  be 
expressed  in  terms  of  the  partial  pressures  of  the  products  and  reactants. 
This  yields 


p,,  = 

a 

^2 

1 CN 

p 

2 

= Of.,  p ... 

^2^2 

2 CN 

P 

2 

= 0*^  P , 

^4^2 

3 CN 

p 

C^N 

= «4  PCN 

(5) 


The  ct  coefficients  are  defined  in  terms  of  the  pressure  equilibrium 
constants  for  the  reactions  (presented  in  Appendix  B)  and  the  vapor  pres- 
sure of  monatomic  carbon.  Therefore,  the  Of  coefficients  are  known  func- 
tions of  the  surface  temperature. 
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The  elemental  mass  fractions  at  the  wall  are  related  to  their  boundary 
layer  edge  values  and  the  surface  mass  transfer  coefficient  by  the  analogy 
between  mass  and  heat  transfer  (assuming  equal  diffusion  coefficients  and 
Lewis  and  Prandtl  numbers  of  unity).  This  yields 


O 1+B' 
w 


1+B' 


"1+B' 


B'  +K, 


P u C-_. 
e e H 


The  partial  pressure  of  any  constituent  is  given  as 


PMK. 


P.  = 
1 


The  set  of  Eqs.  (3)  may  be  expressed  in  terms  of  partial  pressures  by 
multiplying  through  by  the  product  of  pressure  and  molecular  weight.  After 
this  multiplication,  the  equations  are  combined  with  Eqs.  (4),  (5),  and  (7)  to 
yield  a quadratic  for  the  partial  pressure  of  CN.  The  algebraic  details  are 
presented  in  Appendix  C. 


(1  + 28X)(a^  i-a^  +0*3)  (I  + 14X)(l+ct^)  + 

- fp  -LPr-  ) = » 


Xq  ‘'i 

— :z^ — + — ^ — 
16 

e e 


After  the  quadratic  is  solved  for  P^-p^  (use  plus  sign  in  front  of  radical 
for  the  correct  root),  the  remaining  partial  pressures  are  calculated  from 
Eq.  (5).  The  mass  transfer  parameter  is  obtained  from  the  elemental  car- 
bon mass  balance  of  Eq.  (3)  and  is  obtained  as  (see  Appendix  C) 


B'  = B -K^  + 

e 


\ [(12+24a4)Pc^  + (24a2  + ^«“3^^CN  ^E^C.^C.] 

14(l+a4)P^j^  + 28(a^+a2+a3)p2^ 


K 


If  it  is  desired  to  examine  the  relative  contribution  of  the  carbon- 
nitrogen  species,  the  appropriate  cx  coefficients  of  Eq.  (5)  may  be  zeroed. 

If  all  carbon-nitrogen  reactions  are  deleted,  the  system  of  equations  is  solved 
for  the  partial  pressure 

p-EPc. 

Pn,  ' 1 + 28X 


Thus,  for  no  carbon-nitrogen  reactions,  the  mass  transfer  parameter  is 


B'  = B - K + 
e 


(1  +28X) 

e ^i  ^i 

28(p-EPc,) 


Note  that  Eq.  (12)  is  the  analog  of  Bartlett's  only  model  for  an  arbitrary 
number  of  sublimation  species. 

If  there  is  no  nitrogen  in  the  atmosphere  = oj  , Eq.  (12)  yields 


B'  = b''  - K + ( ^ + 1^)  

^1/  P 


After  B'  is  solved  for,  the  molecular  weight  at  the  wall  is  calculated 
from  the  carbon  mass  balance  of  Eq.  (3). 

(b'  -B%K^  ) = 12(1+2c.,)P^^  + 24(c.2+^<^>PcN  + 2:Mc.Pc. 


After  the  molecular  weight  from  Eq.  (14)  is  obtained,  the  species  mass 
fractions  may  be  calculated  from  the  relationship  between  partial  pressure 
and  mass  fraction  [Eq,  (8)]. 

A closing  comment  is  in  order.  Notice  that  the  solution  could  be 
generalized  by  including  the  reduction  of  CO^  to  CO.  This  would  allow  the 
calculation  of  diffusion  controlled  recession  at  low  temperatures.  However, 
at  low  temperatures,  the  recession  is  likely  to  be  rate  limited,  and  the 
diffusion  controlled  calculated  would  provide  only  an  upper  bound  for  the 


w 


actual  mass  loss.  Furthermore,  the  inclusion  of  the  reduction  reaction 

would  yield  a quartic  to  be  solved  for  either  the  CN  or  CO  partial  pressure. 
Although  an  analytic  solution  for  a quartic  is  available,  it  is  considerably 
more  complex  than  evaluating  a quadratic.  For  these  reasons,  the  CO^ 
reaction  was  not  included. 

As  mentioned  earlier,  the  results  of  this  section  may  be  used  for 
carbon  sublimation  in  air.  For  Argon  free  air,  the  following  values  are  used: 


Kq  = 0.232 
e 


=0.768 

e 


Kc  =0 
e 


Kl  = 0 

e 


B = 0. 174 


X = 0. 01888 


9’/-£cediA'i^  ~ 
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SECTION  III 

CARBON  SUBLIMATION  IN  A C-H-O-N-I  ATMOSPHERE 
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In  this  section,  the  sublimation  of  carbon  is  examined  in  an  atmosphere 
whose  elemental  composition  consists  of  carbon,  hydrogen,  oxygen,  nitrogen 
and  an  inert  monatomic  species.  The  applications  of  interest  for  this  atmos- 
phere consist  of  planetary  entry  cases  (e.g.,  Jupiter)  and  ablation  testing  in 
rocket  engines.  Since  the  planetary  entry  application  involves  much  higher 
surface  temperatures  than  the  rocket  engine  environment,  the  analysis  is 
divided  into  high  and  low  temperature  domains. 

The  inclusion  of  hydrogen  in  the  elemental  composition  is  accommodated 
in  the  formulation  by  ignoring  the  carbon-nitrogen  reactions.  That  is,  since 
carbon  reacts  preferentially  with  hydrogen,  the  carbon-nitrogen  reactions  are 
assumed  to  exert  only  a second-order  influence  on  the  surface  mass  transfer 
coefficient. 

The  surface  temperature  and  pressure  and  the  boundary  layer  edge 
elemental  composition  are  presumed  known.  As  before,  the  partial  pres- 
sures of  the  carbon  sublimation  species  are  equated  to  their  corresponding 
vapor  pressures  and  are  then  known  functions  of  temperature.  Again,  the 
elemental  mass  fractions  at  the  boundary  layer  edge  must  sum  to  unity 

+K^  + Kq  + Kj  =1  (16) 

e e e e e 

The  high  and  low  temperature  regimes  are  considered  below. 

A.  HIGH  TEMPERATURE  REGIME 

The  surface  temperature  is  presumed  to  be  high  enough  so  that  all 
oxygen  is  bound  in  CO.  Therefore,  CO^  and  H^O  exist  only  in  trace  amounts 
and  are  ignored.  Although  the  surface  temperature  is  such  as  to  preclude  the 
formation  of  water  vapor,  it  is  assumed  to  be  insufficient  to  dissociate  N^  or 
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CO.  Thus,  monatomic  nitrogen  and  oxygen  are  not  present  at  the  wall.  The 
hydrogen  dissociation  reaction  is  included,  along  with  the  two  dominant 
carbon-hydrogen  reactions  and  an  arbitrary  number  of  carbon  sublimation 


species . 


The  gas  phase  reactions  are  specified  as 


C2H-^H  + 2C 


C^H^-*2H  + 2C 


The  elemental  mass  balances  and  molecular  weight  at  the  wall  are  given  by 


K - JAk 

~ 28  CO 

" 28  ^CO  2^  25  ^ 


Ki  = Kj 

w 


± = ^CO  /^^2  ^»2 

M 28  28  ^ 2 ^ ^ 


25 


_L  +y^ i 


The  relationship  between  the  wall  elemental  mass  fractions  and  the 
boundary  layer  edge  values  is  supplied  by  Eq.  (7),  with  the  exception  of 
hydrogen,  which  is  given  as 


H 1 +B' 
w 
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The  partial  pressures  for  the  reactants  and  products  of  Eq.  (16)  are 
given  by 


C^H 


- 


> s a P 

C2H2 


2 

a 

7 H 


(20) 


The  01  coefficients  are  defined  in  terms  of  the  equilibrium  constants  and 
the  vapor  pressure  for  monatomic  carbon. 


(21) 


Equations  (18)  are  multiplied  by  the  product  of  pressure  and  molecular 
weight  to  cast  them  in  terms  of  partial  pressures.  Then,  combining  Eqs.  (7), 
(18),  (19),  and  (20)  yields  a quadratic  for  the  partial  pressure  The 

algebraic  operations  are  similar  to  those  of  Appendix  C. 


^2<“6"“7'Ph  + ^I<*"“5'Ph  -(P-SPc.)^H  = ° 


K 


P - 


N 


K 


28 


O 

e 

TF 


K, 


M, 


e 

^2  = 2/3  + (22) 
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After  solving  the  quadratic  for  all  the  remaining  partial  pressures 
are  easily  evaluated.  The  surface  mass  transfer  parameter  is  evaluated 
from  the  carbon  mass  balance  of  Eq.  (18)  and  yields 


B'  = B - K + ^ 

e P.£p^_.(«^+a^)P^-(l+«3)P^ 


(23) 


As  mentioned  earlier,  the  validity  of  this  expression  for  B'  presumes 
that  the  carbon-nitrogen  reactions  [which  were  ignored  in  obtaining  Eq.  (23)] 
are  insignificant  in  comparison  with  the  carbon-hydrogen  reactions.  In  the 
limit  that  the  elemental  mass  fraction  of  hydrogen  ^ Kj^  ^ vanishes,  this 
assumption  will  become  progressively  worse.  However,  if  only  trace  amounts 
of  hydrogen  were  present,  it  could  be  ignored  and  the  mass  transfer  coefficient 
evaluated  from  the  results  of  Section  II. 

In  any  event,  notice  that  the  carbon-nitrogen  reactions  were  deleted 
for  simplicity  and  not  from  necessity.  That  is,  it  is  possible  to  include  all 
the  carbon-nitrogen  reactions  of  Eq.  (2)  along  with  the  carbon-hydrogen 
reactions  of  Eq.  (17)  and  still  obtain  an  analytic  solution  for  the  mass  trans- 
fer coefficient.  However,  the  determination  of  the  constituent  partial  pres- 
sures would  require  the  solution  of  a quartic.  Although  a quartic  does  possess 
an  analytic  solution,  it  was  felt  that  the  additional  complexity  was  not  war- 
ranted for  the  intended  application. 

The  molecular  weight  is  obtained  from  either  the  hydrogen  or  carbon 
mass  balances  [Eq.  (18)].  The  hydrogen  mass  balance  yields 

1 +B'  ^H  ~ 
e 
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In  the  previous  subsection,  the  surface  temp)erature  was  assumed  to  be 
high  enough  so  that  the  water  gas  reaction  was  insignificant.  At  lower 
temperatures,  however,  the  formation  of  water  vapor  and  carbon  dioxide 
has  a first-order  effect  on  the  surface  mass  transfer.  Since  this  regime 
may  be  important  for  ablation  testing  in  rocket  engines,  it  is  considered  here. 

The  surface  temperature  is  assumed  low  enough  so  that  atomic  oxygen, 
hydrogen,  and  nitrogen  are  not  present  at  the  wall.  All  the  oxygen  is 
assumed  to  be  bound  up  in  H^O,  CO^  and  CO.  Again,  an  arbitrary  number  of 
carbon  sublimation  species  are  included,  as  well  as  a single  carbon-hydrogen 
reaction.  The  reactions  are 


H^O  + CO  -►CO^  + 


C + CO^-^Z  CO 


+2C 


(25) 


The  equilibrium  relationships  relate  the  partial  pressures  of  products 
and  reactants. 


P = a p p 

^H^O  8 CO 


P = a P 
CO^  9 CO 


p = a p 

C^H^  lO^H^ 


(26) 
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The  01  coefficients  are  given  by 


i 


I 

] 

I 


I 

i 


Kp  Pc 

^8  9 


a 


10  K, 


10 


The  mass  balances  and  molecular  weight  are  given  by 


K +11k  +iiK 

“ 28  ^CO  44  co^  18  H^O 


K = K 

J-W  ^ 


w 


’ 18  H^O 


“ 28  ^CO  44  ^CO-  2^ 
w 2 2 2 1 


J_ 

M 


K 


K 


CO 


N, 


K. 


K 


C^H 


28 


CO  O K 

2E"  + nx“  + 


2 2 


K 


C. 


M 


C. 


(27) 


(28) 


Equations  (7) 
and  boundary  layer 


and  (19)  relate  the  elemental  mass  fractions  at  the  wall 
edge. 
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If  Eqs.  (28)  are  cast  into  partial  pressures,  the  system  may  be  re- 
arranged to  yield  a cubic  for  the  partial  pressure  of  CO. 

APcO  " BPco  + CPco  + D = » 


^ + 4Xj  +ag(l +a^Q)^112K^  + 7Kq  + 8X 

C = (1 +a^Q)/56Kj^  + 7Kq  t 4x\  -t 


+ 7a8f8K^^-Ko  JfP-EPc 


D = - 7K_  P - 


X = K-,  + ■ .. 

N 

e I 


After  the  cubic  is  solved,  the  partial  pressure  of  is  obtained  from 


28K_  P-Ep 


I6X  + 28K0  iP^o 

< e/ 


32X+28Ko^^gP^O 


28Kq  +ag/l6X  + 28KQ\p^O  + 28ajoKQ 


The  carbon  mass  balance  is  solved  for  B'  to  yield 


B'  = -Kc  4 
e 


® j^^CO  "^^9^00  I2I 

^CO  ^°‘9^CO  “s^CO^H. 


i i- 


i 


1 

I 


It  should  be  emphasized  that  the  result  of  Eq.  (31)  was  obtained  from 
the  assumption  that  the  surface  recession  is  diffusion-controlled.  In  actual 
fact,  at  sufficiently  low  temperatures,  the  recession  may  be  rate- limited. 
Therefore,  Eq.  (31)  should  be  viewed  as  providing  an  upper  bound  for  the 
recession  at  low  temperatures. 

After  the  calculation  of  B',  the  molecular  weight  may  be  calculated 
from  one  of  the  mass  balances  of  Eq.  (28).  The  hydrogen  mass  balance 
yie  Ids 

' ^'‘+“iO+“8Pco'Ph, 

e 2 


;i 

1 

,1 
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SECTION  IV 


r 


j 

,i 

:i  COMPARISON  WITH  EXACT  SOLUTIONS 

' In  this  section,  the  analytic  solutions  of  Sections  II  and  III  are  compared 

with  exact  results  obtained  from  the  ACE  (Ref,  2)  thermochemistry  computer 
code.  This  code  solves  the  general  heterogeneous  system,  for  a large  num- 
ber of  species,  with  an  iterative  numerical  scheme  and  is  generally  regarded 
as  the  standard  of  comparison  in  the  reentry  community. 

Ideally,  the  comparisons  between  the  ACE  results  and  the  present  work 
should  be  carried  out  using  identical  assumptions  and  thermochemical  data. 
Thus,  the  ACE  results  were  generated  using  the  option  of  equal  diffusion 
coefficients.  Two  sources  of  thermochemical  data  [Horton  (Ref.  3)  and 
JANNAF  (Ref.  4)]  were  used  to  generate  the  analytical  solutions,  whereas 
only  JANNAF  data  are  used  in  ACE.  However,  care  was  taken  to  insure  that 
the  data  for  the  major  species  were  identical  to  JANNAF.  The  slight  dif- 
ferences among  the  secondary  species  are  not  significant  enough  to  exert 
any  impact  on  the  mixture  properties  or  mass  transfer  coefficient. 

The  first  set  of  comparisons  is  for  carbon  sublimation  in  air  and  for 
comparison  of  the  analytic  solutions  of  Section  II  with  ACE  results.  The 
values  of  Eq.  (15)  were  used  to  represent  air,  and  the  first  three  carbon  sub- 
limation species  were  included  in  the  analytic  solution.  Figures  1 and  2 show 
the  mass  loss  coefficient  (B')  as  a function  of  surface  temperature  for  two 
pressures.  The  two  solutions  of  Section  II  are  shown  as  solid  lines  and 
indicate  the  influence  of  the  carbon-nitrogen  reactions  on  the  mass  transfer. 
The  upper  curve  [Eq.  (10)]  includes  all  the  dominant  carbon-nitrogen  reactions 

3 

Horton,  T.E.,  and  W.A.  Menard,  "A  Program  for  Computing  Shock-Tube 
Gas  Dynamic  Properties,  " Jet  Propulsion  Laboratory,  Pasadena,  California, 
TR-32-1350,  January  1969. 

4 

JANNAF  Thermochemical  Tables,  National  Bureau  of  Standards,  NBS-37, 
June  1971. 
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and  is  in  excellent  agreement  with  the  ACE  results.  The  lower  curve  [Eq. 
(12)]  ignores  carbon-nitrogen  reactions  and  is  in  considerable  error.  It 
seems  apparent,  therefore,  that  ignoring  carbon-nitrogen  reactions  (at  least 
for  the  1969  JANNAF  data)  is  a poor  assumption.  The  distribution  of  species 
at  the  wall  is  shown  in  Figs.  3 and  4.  The  solid  lines  are  obtained  from 
Eqs.  (5),  (8),  and  (9).  Again,  the  agreement  with  ACE  is  excellent  for  the 
major  species  and  good  for  the  secondary  species.  As  mentioned  above, 
the  slight  discrepancy  between  the  secondary  species  is  due  to  differing  sets 
of  thermochemical  data. 

The  next  set  of  comparisons  is  for  a carbon  surface  ablating  in  a rocket 
engine  exhaust.  In  particular,  the  RPL  (Air  Force  Rocket  Propulsion 
Laboratory)  rocket  engine  is  used  as  the  test  case.  The  engine  burns  benzo> 
itrile  and  LOX,  at  a mixture  ratio  of  2.  8,  and  yields  the  following  elemental 
composition: 


e 

= 0. 012755 

e 

= 0.  03571 

^c 

e 

= 0. 2143 

e 

= 0.  7372 

^i 

= 0 

(33) 


The  surface  mass  transfer  parameter  from  Section  III  is  shown  com- 
pared to  the  ACE  results  for  two  pressures  in  Figs.  5 and  6.  Both  the  high 
and  low  temperature  expressions  [Eqs.  (23)  and  (31)]  are  plotted,  and  both 
are  in  excellent  agreement  with  ACE  values  in  their  respective  domains  of 
validity.  That  is,  the  low  temperature  result  of  Eq.  (31)  (shown  as  a solid 
line)  agrees  well  with  the  ACE  results  for  low  to  moderate  temperatures,  and 
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T 


the  high  temperature  result  of  Eq.  (23)  (shown  as  a dashed  line)  is  in  good 
agreement  for  moderate  to  high  temperatures.  It  is  believed  that  the  excel- 
lent agreement  shown  in  Figs.  5 and  6 justifies  the  assumption  that  carbon- 
nitrogen  reactions  may  be  ignored  relative  to  the  carbon-hydrogen  reactions. 
Notice  that  the  low  temperature  result  has  the  correct  behavior  at  high 
temperatures,  due  to  the  inclusion  of  the  carbon  sublimation  species,  but 
does  not  rise  fast  enough  because  C2H  was  ignored  in  the  interests  of  ob- 
taining an  analytic  solution. 

The  ACE  results,  in  keeping  with  the  assumptions  of  this  report,  are 
all  for  diffusion-controlled  surface  reactions.  Thus,  at  low  temperatures, 
where  the  reactions  maybe  rate  limited,  the  results  of  Figs.  5 and  6 should 
be  interpreted  as  upper  bounds  to  the  actual  recession  rate. 
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SECTION  V 

SUMMARY  AND  CONCLUSIONS 


A new  class  of  analytic  solutions  has  been  presented  for  mass  transfer 
from  a carbon  surface  subliming  into  atmospheres  of  varying  elemental 
composition.  The  solutions  are  computationally  trivial,  requiring  only  the 
evaluation  of  a single  quadratic  equation  (in  the  diffusion  limited  regime)  and 
are  in  excellent  agreement  with  exact  results  generated  by  the  ACE  code. 

For  an  atmosphere  of  elemental  C-N-O-I,  a single  analytic  solution  for 
B'  was  obtained  in  the  temperature  regime  where  the  oxidation  reaction  is 
diffusion-limited.  An  arbitrary  number  of  carbon  sublimation  species  is 
included  in  the  solution,  as  well  as  all  the  dominant  carbon-nitrogen  reactions. 

When  hydrogen  is  added  to  the  elemental  composition,  the  analytic 
solution  is  achieved  by  ignoring  carbon-nitrogen  reactions  and  treating  the 
dominant  carbon-hydrogen  reactions.  Two  solutions  are  presented  for  this 
case:  one  valid  in  the  temperature  domain  where  the  oxidation  reaction  is 
diffusion-limited  and  a second  valid  at  lower  temperatures  where  the  water 
gas  reaction  is  dominant. 

The  solutions  should  be  particularly  suited  for  applications  in  which 
the  atmospheric  elemental  composition  is  frequently  changed  and  in  situations 
where  a knowledge  of  the  analytic  behavior  of  B'  with  temperature  is  desired. 
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APPENDIX  A 

VAPOR  PRESSURE  OF  CARBON  SUBLIMATION  SPECIES 


As  discussed  in  the  main  text  (Section  II),  the  vapor  pressure  of  each 
gaseous  carbon  species  is  assumed  to  be  a known  function  of  temperature. 
These  vapor  pressures  are  presented  in  this  appendix. 

The  gaseous  carbon  species  in  equilibrium  with  the  solid  carbon  sur- 
face are  denoted  by  C^^.  The  vaporization  reaction  is  given  as 


iC  ►C. 

solid  1 


For  this  reaction,  the  third  law  gives 


(A-1) 


Ixv  P, 


fO-hJ 


(A-2) 


solid 


where  Ert(C.)  is  the  heat  of  formation  of  C.  at  0 K and  the  heat  of  formation  of 
0 1 1 

solid  carbon  has  been  set  equal  to  zero  by  convention. 

The  free-energy  function  and  partition  function  for  the  gaseous  species 
are  related  as  follows: 


r»-HS 


- in  Qp(C.) 


(A-3) 


The  free-energy  function  is  related  to  the  specific  heat  by 


fO-h? 


-dT  - s: 


(A-4) 
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Note  that  the  JANNAF  (Ref.  4)  thermochemical  data  uses  a base 
temperature  of  298°K.  The  free-energy  functions  at  0°K  and  298°K  are 
’ Uced  by 


0"  298 


(A-5) 


The  free-energy  function  referenced  to  298°K  and  the  difference  in 
enthalpies  between  0°K  and  298°K  are  both  tabulated  in  JANNAF. 

For  solid  carbon,  the  JANNAF  data  yields 


F -Hq  ^ F -H^gg  252 


(A -6) 


The  curve  fits  for  specific  heat  in  the  ACE  (Ref.  2)  code  were  used  to 
evaluate  the  free-energy  function  for  solid  carbon  by  using  Eq.  (A-4).  The 


final  result  is 


= A(l-/nT)  - y BT  - - § + E 

^ 2T^ 


(A-7) 


T < 3000°K: 


5.  861 

0.  954  X 10“'^ 
-0. 7666  X 10^ 
4323. 3 
35. 235 
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T > 3000°K: 

A = 4. 85 

B = 0.  2916  X 10"^ 

C = 0.  3072  X lO"^ 

D = 893. 7 
E = 27.516 

where  all  temperatures  are  in  degrees  Kelvin  and  the  units  of  the  free-energy 
functions  are  cal/mol-°K. 

The  use  of  Eqs.  (A-6)  and  (A-7)  in  Eq.  (A-2)  provides  all  the  necessary- 
information  about  the  solid  phase. 

The  partition  functions  were  used  to  obtain  the  free-energy  function  for 
C and  and  are  presented  in  Appendix  B.  The  free-energy  function  was 
again  obtained  from  the  specific  heat  curve  fit  in  the  ACE  code  (196  9 JANNAF). 
The  free-energy  function  is  identical  to  the  form  of  Eq.  (A-7)  with  the 
following  values  for  the  constants: 

HQ-H298  = - 2811  cal/mol 

T < 2500°K: 

A 
B 

i C 

D 
E 

T > 2500°K: 

A 
B 
C 
D 
E 


= 10.223 

= 1. 0066  X 10’^ 

= -0. 44592  X 10^ 
= 4587.  525 
= 3. 722 

= 12. 524 
= 2.  3055  X lO""^ 

= -0.  27025  X lo"^ 
= 8817. 501 
= 19.966 
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APPENDIX  B 

EQUIUBRIUM  CONSTANTS 


This  appendix  contains  the  necessary  information  required  to  calculate 
the  pressure  equilibrium  constants  which  are  used  in  the  main  text. 

The  reactions  considered  are  all  of  the  following  form: 


(B-1) 


where  the  A^,  are  the  reactants  and  the  B^  the  products.  The  equilibrium 
constant  for  this  reaction  is  defined  in  terms  of  the  partial  pressures  as 


^1  ^2 
P P ^ 

®1  ®2 

Kp  = 

^ ^1  ^2 
Pa  Pa 

Ai  A^ 


(B-2) 


where  P.  denotes  the  partial  pressure  of  A.. 

Af 

The  equilibrium  constant  may  be  computed  from  a knowledge  of  the 
partition  function  (or  equivalently,  the  free  energy  function)  of  each  species 
involved  in  the  reaction.  That  is 


log  Kp  = + b^log  Q (B^  + b^log  Q (B^) 


ailog  Q (Aj)  - a^log  Q (A^) 


where  E^  is  the  heat  of  formation  at  0 K. 
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The  equilibrium  constants  are  defined  for  the  indicated  reactions  as 
follows: 


(N^ 

+ 2C-*>2  CN) 

(^2 

N2-^2  CN) 

(C4 

N2->2  CN  + 2C) 

(^2 

N-^CN  + C) 

(C2 

H-*-H  + 2C) 

(C2 

H2-^2H  + 2C) 

■^7 

(H2 

-►2H) 

H,( 

D + CO  -►CO^  + 

9 

(C  • 

f CO^-^Z  CO) 

(C 

+ 2C) 

The  heat  of  reaction  term  is  defined  for  each  reaction  as  follows: 


i 

i 


1 


K, 


4Eq/R 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


-67260.  0 
66993.  0 

211521. 0 

71143. 0 
138419.  0 
195816.  0 

51984. 0 
-4861.  0 

-65687.  0 
143832.  0 


K 


The  partition  functions  used  in  this  report  were  obtained  from  the  data 
of  Horton  (Ref.  3)  and  JANNAF  (Ref.  4)  and  are  given  in  Table  B-1.  The 
partition  function  may  be  calculated  from  the  free  energy  functions  presented 
in  Appendix  A.  Similarly,  the  free  energy  function  for  is  also  obtained 

from  Eqs.  (A-5)  and  (A-7)  of  Appendix  A,  with  the  following  constants. 


T < 2500°K: 


Hq-H298  = - 4239  cal/mol 


A = 28. 903 
B = 0.  13659  X 10-2 
C = -0.  14872  X 10'7 
D = 1.  3664  X lO'^ 

E = 102.664 


T > 2500°K: 


A 

= 32 

. 053 

B 

= 0. 

109  X 10- 

3 

C 

= -0 

. 15347 

X 

10 

D 

= 1. 

763  1 X 

10 

4 

E 

= 124. 171 

48 


Table  B-1.  Partition  Functions 


log  Q (CN)  = 3.5  log  T + 1.  147  + log  2.0/[2.  72  - 2.  72  exp(-2958/T)]  + 4 exp(- 1 3 1 14/T) /[2.  46  - 2.46  exp 


APPENDIX  C 

DETAILS  OF  THE  ANALYTIC  SOLUTION 


The  algebraic  operations  leading  to  Eqs.  (9)  and  (10)  of  Section  II  are 
presented  in  this  appendix. 


As  indicated  in  the  text,  Eq.  (3)  is  cast  in  terms  of  partial  pressures 
by  multiplying  through  by  the  product  of  pressure  and  molecular  weight. 
This  operation,  in  conjunction  with  the  use  of  Eq.  (7),  yields  the  following: 


PM  ~ - 16,  p 

1 +B'  ^O  ^CO 

e 


TW  «N  = Pn,  + PcN  * Pc,N,  * '’C,N 

e 2 2 2 4 2 2 


12  PCO  + ‘2  PcN  * 24  + 48  P 

2 2 4 2 


+ 24Pc^N+  LMc.Pc, 


1 1 


PM 

iw^i  =^1^1 

e 


(C-l) 


The  sum  of  the  partial  pressures  must  equal  the  mixture  pressure. 
That  is 


P * P^,..  + P..  + P^„  + P^  + P^  + P^  .,  + Pj  + £Pj3  (C-2) 


CO  ^ N^  ■ ^ CN  - C^N^  C4N2  C^N 
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Substitution  of  P^q  and  from  Eq.  (C-1)  into  Eq.  (C-2)  yields 


MP 

r+B^ 


e_  e 

TT^  Mj 


■ ^ ‘^C.  ■ ■ ^CN 

1 2 


- P - P - P 

C,N-  C.N-  C,N 

2 2 4 2 2 


(C-3) 


Combining  Eq.  (C-3)  with  the  nitrogen  mass  balance  of  Eq.  (C-1) 


yields 


(28X+1) 


2 ^ ^^2^2  ^ ^ 


+ (14X+ 1)  + = P 


{C.4) 


Substitution  of  the  equilibrium  relationships  of  Eq.  (5)  into  Eq.  (C-4) 
yields  the  result  of  Eq.  (9). 

Substitution  of  iJ^to  the  carbon  mass  balance  of  Eq.  (C-1)  yields 


-PM.  / 
1 +B' 


B'  + - B = 12  P^,,  + 24  P^  + 48  P^  + 


CN 


^C2N  P^ 


(C-5) 


1 1 


The  term  PM/ (1  +B')  may  be  eliminated  from  Eq.  (C-5)  by  using 
either  Eq.  (C-3)  or  the  nitrogen  mass  balance  of  Eq.  (C-1).  Equation  (10) 
was  obtained  by  using  the  latter  and  the  equilibrium  expressions  of  Eq.  (5). 
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APPENDIX  D 


GENERALIZED  ANALYSIS  OF  THE  C-N-O-I  SYSTEM 


I 


i-i 


The  analysis  of  Section  II  ignored  the  presence  of  CO2  at  the  wall  and 
as  sunned  that  all  the  oxygen  was  bound  in  CO.  This  limits  the  validity  of 
the  results  to  surface  temperatures  high  enough  to  preclude  the  formation  of 
CO^.  As  mentioned  on  page  14,  at  the  low  surface  temperatures  required 
for  the  existence  of  CO^  the  surface  recession  is  likely  to  be  rate  limited, 
and  the  diffusion-controlled  calculation  provides  only  an  upper  bound  to  the 
actual  recession.  There  are,  however,  situations  in  which  a knowledge  of 
the  upper  bound  recession  suffices,  and  the  complexity  of  the  surface  kinetic 
reactions  may  be  ignored.  In  order  to  obtain  a reasonable  upper  bound 
value  for  B',  the  CO^  reaction  must  be  included  in  the  formulation.  It  is  the 
purpose  of  this  appendix  to  generalize  the  analysis  of  Section  II  by  including 
CO^  in  the  formulation. 

The  analysis  is  considered  below  in  two  parts.  In  the  first  part,  the 
problem  is  formulated  and  discussed.  Since  the  generalized  formulation 
leads  to  a quartic  equation  (which  is  time  consuming  to  evaluate),  an  approxi- 
mate (but  very  accurate)  method  of  solution  is  presented  in  part  2. 

1.  FORMULATION 

The  reactions  include  those  of  Eq.  (2)  plus  the  CO2  reduction  reaction 
which  is  given  as 

C + CO^-^Z  CO  (D-1) 

The  equilibrium  relationships  include  those  of  Eq.  (5)  as  well  as  the 
relation  between  the  CO  and  CO^  partial  pressures. 


P = a P 
CO2  5 ^CO 


“5  “ P^K 


(D-2) 


C"P, 


Note  that  the  coefficient  of  Eq.  (D-2)  is  not  the  same  as  that  of  Eq.  (21). 
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The  elemental  mass  balances  of  Eq.  (3)  have  the  additional  terms 
added  as  follows; 


K - K + K 
0,„  ■ 28  CO  44  CO- 

W L 


^ 12  12  12  24 

K_  + 

C,,,  28  CO  44  CO-  26  CN  52  C^N- 

W c c i. 


H 38  ^C  N "^S^C. 

4 2 2 1 


(D-3) 


The  nitrogen  and  inert  elemental  mass  balances  of  Eq.  (3)  are  unchanged. 
The  molecular  weight  Eq.  (4)  becomes 


1 _ *^CO 
M “ 28 


^CO 

^ 1 

“"Na  . ^cN 

44  ^ 

28  26 

^C  N 
^4  2 

^C,N  K 

1 ^ \ ^ 

76 

38  M, 

C^N 

+ + 


(D-4) 


Equations  (7),  relating  the  wall  and  boundary  layer  edge  elemental 
mass  fractions,  are  unchanged. 

In  the  same  manner  as  in  Section  II,  the  set  of  equations  may  be 
manipulated  to  yield  the  following  equations  for  the  partial  pressures  of  CO 
and  CN. 

^“1  ■'^“2  '*’®3^^CN  '*■  2 ■^“4^^CN  ~ ^^\^CO  ^^^l^S^CO 


MfKo 


(D-5) 
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- (“l  +“2+“3'PcN-  " +“4>PcN 


Equations  (D-5)  and  (D-6)  may  be  combined  to  yield  a quartic  for  the 
partial  pressure  of  CO.  After  the  quartic  is  solved,  the  CN  partial  pres- 
sure is  calculated  from  Eq.  (D-5)  and  the  remaining  partial  pressures  are 
obtained  from  Eqs.  (5)  and  (D-2).  Following  the  determination  of  the  partial 
pressures,  the  molecular  weight  and  B’  are  calculated  from  Eqs.  (D-3) 
(recast  in  terms  of  partial  pressures).  These  equations  yield 


PM 

1+B' 


^CO  " Pco, 


PM 

1+B' 


(b.  4 


Pco  " Pc02  ^CN  + '’C2N  * 


+ PC2N2  " '’c^N2  P 


c. 

1 1 


(D-7) 


This  completes  the  generalization  of  the  analysis  of  Section  II  to 
include  the  dominant  low  temperature  reaction  [i.e.,  Eq.  (D-1)]. 

A comment  is  in  order  concerning  the  advisability  of  including 
among  the  low  temperature  reactions.  That  is,  it  is  quite  easy  to  add  an 
additional  reaction  specifying  equilibrium  between  C,  CO  and  (or, 
equivalently,  O2,  and  CO)  and  including  O^  in  the  oxygen  elemental 

balance  of  Eq.  (D-3)  and  the  molecular  weight  equation,  (D-4).  Inclusion  of 
O^  in  the  analysis  changes  only  the  coefficients  of  the  quadratic  terms 
involving  the  CO  partial  pressure  of  Eqs.  (D-5)  and  (D-6).  Calculations 
were  done  with  the  O^  included,  and  it  was  found  to  exert  no  influence  on  the 
results  down  to  temperatures  as  low  as  600°R.  Therefore,  it  was  deleted 
in  the  above  description  of  the  analysis. 
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The  limiting  behavior  of  the  mass  transfer  parameter  B'  for  large 
and  small  temperatures  is  considered  next.  Since  was  ignored  in  the 
analysis,  as  the  temperature  vanishes,  only  CO^  is  present  at  the  wall  and 
B'  approaches  a constant  value.  That  is 

T— ►O  (CO^  plateau): 


B-  -K 


B'  = 0. 087  for  air 


(D-8) 


where  the  values  of  Eq.  (15)  were  used  for  air.  In  contrast,  the  low  tempera- 
ture limit  for  the  analysis  of  Section  II  (where  only  CO  is  present  as  the 
temperature  vanishes)  is  given  as 

T— ►O  (CO  plateau); 

12  ~ ~ 

=T^^o  -^c 

e e 


B'  = 0.174  for  air 


(D-9) 


There  exists  an  upper  limit  temperature  beyond  which  no  equilibrium 
solution  exists.  This  limiting  temperature  is  a function  of  pressure  and  is 
defined  by  the  sum  of  the  carbon  vapor  pressures  equal  to  the  mixture  pres- 
sure. This  yields 


S^C.  ^^limit^  " 


(D-10) 


r » 


The  high  temperature  behavior  of  the  generalized  system  is  identical  to 
that  of  Section  II  and  is  given  as 

T — -►T  * 

limit 


CN  " (1  + 14A){1  +0^) 
14  B*(p 


CO  12(1+14X)K 


N 


(1  + 


B*  = 


1 1 


00 


(D-ll) 


where  X is  defined  by  Eq.  (9). 


2.  SOLUTION 

As  indicated  above,  Eqs.  (D-5)  and  (D-6)  may  be  combined  to  yield  a 
quartic  for  the  CO  (or  CN)  partial  pressure.  Although  an  analytic  solution 
exists  for  a quartic,  it  is  time  consuming  to  evaluate,  and  the  possibility  of 
more  than  one  plausible  root  complicates  the  issue.  In  order  to  avoid  this 
complexity,  the  solution  for  the  CO  partial  pressure  is  split  into  high-  and 
low-tempe rature  regimes  requiring  only  the  solution  of  a quadratic  in  each 
regime. 

In  the  low  temperature  regime,  the  carbon-nitrogen  species  are  ignored 
and  the  following  quadratic  is  obtained. 


a^Cl  + 32Xj  + 32X2)P^q  + (1  + l6Xj  + IbX^)?^-^^  + 
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There  is  only  one  positive  root,  and  it  is  obtained  by  using  the  positive 
radical  in  the  quadratic  solution. 

The  high  temperature  solution  for  is  obtained  from  the  results  of 

Section  II.  That  is,  after  solving  Eq.  (9)  for  CO  partial  pressure 

is  calculated  from 

r 2 1 

Pco  =-;7^h'“l+“2+"3'PcN+“<‘+“4>PcNj 

e 

The  correct  value  for  P(-|qi  for  any  temperature,  will  be  the  smaller 
of  the  two  values  obtained  from  Eqs.  (D-12)  and  (D-13).  After  obtaining  this 
value  for  P(-q»  Eq.  {D-5)  is  solved  as  a quadratic  for  the  CN  partial  pres- 
sure, and  the  remaining  partial  pressures  are  calculated  from  Eqs.  (D-2) 
and  (5).  Then  Eqs.  (D-7)  are  solved  for  B'  and  the  molecular  weight. 

In  this  manner,  an  approximate  solution  is  obtained  for  the  generalized 
formulation  without  the  need  to  evaluate  a quartic  equation.  The  solution  is 
continuous  for  all  temperatures  and  is  quite  accurate. 

The  calculations  of  Section  II  for  an  ablating  carbon  surface  in  air 
were  recalculated  with  the  analysis  of  this  appendix  over  a much  wider  tem- 
perature range.  At  high  temperatures,  the  results  were  identical  to  those 
shown  in  Figs.  1 through  4.  The  low  temperature  results  for  B'  are  com- 
pared in  Fig.  D- 1 with  solutions  from  the  ACE  (Ref.  2)  computer  code. 
Again,  the  agreement  is  excellent. 


Fig.  D-1.  Carbon 
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NOMENCLATURE 

diffusion  limited  mass  transfer  coefficient  [defined  in  Eq.  (10)] 
mass  transfer  parameter  [defined  in  Eq.  (7)] 
pressure  equilibrium  constant 
mass  fraction  of  element  i 
mass  fraction  of  species  i 
molecular  weight  of  gas  mixture 
molecular  weight  of  inert  species 
molecular  weight  of  species  i 

molecular  weight  of  carbon  sublimation  species 
surface  pressure 
partial  pressure  of  species  i 

vapor  pressure  of  carbon  sublimation  species  C. 
surface  temperature 

Subscripts 

e boundary  layer  edge  values 

w wall  (surface)  values 


B 

B' 

K. 

1 

K. 

1 

M 

Mi 

M. 


M 

P 

P. 


c. 

I 


c. 

1 


